Swarming of Self-Propelled Camphor Boats 
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When an ensemble of self-propelled camphor boats move in a one-dimensional channel, they 
exhibit a variety of collective behaviors. Under certain conditions, the boats tend to cluster together 
and move in a relatively tight formation. This type of behavior, referred to as clustering or swarming 
here, is one of three types recently observed in experiment. Similar clustering behavior is also 
reproduced in simulations based on a simple theoretical model. Here we examine this model to 
determine the clustering mechanism and the conditions under which clustering occurs. We also 
propose a method of quantifying the behavior that may be used in future experimental work. 

PACS numbers: 05.65.-|-b,45.50.-j 
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I. INTRODUCTION 

Ensembles of many motile, interacting bodies may ex- 
hibit a variety of organized phenomena. One such behav- 
ior is the grouping of several individuals into a spatially 
compact formation, which is maintained as the bodies 
move. We refer to this type of collective motion as clus- 
tering. It has been observed in the motion of ants [l|, Q 
and numerous other biological systems [^-(^ as well as 
some artificial systems In many cases, the individ- 
ual members of the group do not interact directly, but 
rather change their surroundings in ways that influence 
the behavior of other members. Camphor boats belong 
to this category because, aside from collisions, they in- 
teract only through the background camphor field. An- 
other example of particular relevance to this work is the 
pheromone trail created by an ant, which signals other 
ants to follow. The experimental observations of John 
et al. included moving clusters of ants, which they 
labeled platoons. These platoons share many character- 
istics with the clustering of camphor boats as described 
further below. 

Camphor boats provide a simple way to experimen- 
tally examine collective motion of self-propelled elements. 
Recent camphor boat experiments have displayed vari- 
ous types of collective behavior such as two-boat interac- 
tions 0,1^ and traffic jams |9|]. Two-dimensional systems 
of cam pho r particles have also displayed pattern forma- 
tion [6l,|lC| , but these systems are fundamentally different 
in that they contain free camphor particles with no di- 
rectional orientation, whereas camphor boats are light, 
plastic disks with attached camphor particles that are 
driven in a specified direction by surface tension gradi- 
ents. 

In previous work [9| we presented three types of col- 
lective behavior. Primarily, the boats displayed either 
roughly homogeneous, constant velocity flow or con- 
gested flow(FiglT^ and[T|D) depending on the number of 
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FIG. 1. Space-time diagrams from experimental data show- 
ing the three modes of collective behavior: a)homogeneous(35 
boats, 25° C) b)jammed(40 boats, 25° C) c)clustered(10 boats, 
50° C). Black lines indicate the trajectory of boats. 



boats in the system. The latter occurs with a larger 
number of boats, and is characterized by periodic ac- 
celeration and deceleration. In addition, we observed a 
mode in which the boats gathered together into high den- 
sity groups that were maintained over time as shown in 
FigHJ;. This clustering behavior was briefly described 
therein, but the mechanism and necessary conditions for 
clustering remained under investigation. The theoretical 
analysis presented here reveals the mechanism and nec- 
essary conditions, and provides a method of quantifying 
clustering. In this paper, we will describe this theoret- 
ical mechanism for clustering and present the results of 
numerical simulations. 

We will deflne clustering for this study as the behavior 
in which two or more moving boats gather into a rel- 
atively close formation and remain in that formation as 
they continue moving. This is distinct from the congested 
behavior previously studied [oj in that the boats do not 
experience periodic acceleration, and the formation per- 
sists for a long time. Also, the phenomenon can occur 
with a small total density, whereas congestion typically 
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occurs under higher densities. 



II. MODEL 

There have been several similar models proposed to de- 
scribe the motion of a camphor boat 0, M, Ell, [III • These 
models share the same basic concepts, but differ slightly 
in the mathematical details. The model presented here 
borrows ideas from these. The camphor boats are thin, 
circular, plastic disks with a smaller camphor pellet at- 
tached to the underside so that half of the pellet overlaps 
the edge of the disk and half extends beyond it. We will 
name this overlapping region the back of the boat, and 
the opposite side is the front. All boats are oriented in 
the same direction. If the boats are confined to a long, 
narrow channel in which they cannot pass each other or 
turn around, the motion is essentially one-dimensional. 
The 1-d equations of motion of a boat are given by eqU] 
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where m is the mass of the boat, L is the diameter, and 
fi is the viscosity constant of the water. The position, 
Xj, is defined as the back edge of the boat, which is also 
the center of the camphor pellet. The second term on 
the right represents the difference in surface tension be- 
tween the front and back of the boat as a function of the 
camphor concentration, c{xi + L) and c{xi) respectively. 
7(c) is typically approximated by a sigmoidal function. 
We will use the relation in eql^] 



7(c) 



water ca7nphor 



1 



'^camphor 



(2) 



where ^water = T^g/s^ and ^camphor = 505/5^ are the 
surface tension of pure water and a camphor saturated 
solution respectively. The point c = is the midpoint 

between "f water and -fcaraphor- 

The concentration of camphor on the surface of the 
water is constantly changing due to several processes. 
For a system with N boats, it can be approximated by 
the following reaction-diffusion equation |3l 

N 



~kc + a^F{x~x,) (3) 
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dt dx^ 
F{x) = 1 : for |x| < rp, : otherwise 



Here, D is the diffusion constant, fc is a constant com- 
bining the effects of evaporation and dissolution, and 
aF(x — Xi) represents the addition of camphor by each 
boat's pellet which is centered at the point Xi and has ra- 
dius To- To non-dimensionalize the problem, first define 
the following dimensionless quantities. 
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FIG. 2. a) The difference in surface tension increases with 
increasing average c up to Ct. b) Steady state c(x) for one 
boat(dashed) and two boats(solid) moving at constant veloc- 
ity to the right. The peaks approximate the rear edges of the 
boats, which have length 0.6 



Then the dimensionless parameters of the system become 

.' = ^,fc' = ^,r.^:!^-^ (5) 
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Dropping the ' marks, the non-dimensional equations are 
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dt 
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F{x) = 1 : for |x| < iq, : otherwise 
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III. ANALYSIS 

The mechanism responsible for clustering can be found 
in the surface tension function, 7(c). For concentrations 
less than a critical value, q — I/(\/3/3), the slope of 7(c) 
becomes steeper with increasing c. This means that the 
difference in 7 for a fixed difference in c will increase with 
increasing average c up to ct- This is shown in Figl2^. 
To relate this to a group of boats, if the average cam- 
phor level around the leading boat is lower than around 
trailing boats, the leading boat may feel a weaker driv- 
ing force. The trailing boats would then move at a higher 
velocity and approach the leading boat. 

We can easily apply this idea to the case of two boats. 
First, assume that the boats are moving with constant 
velocity, v, and that the camphor concentration in the 
moving reference frame of the boats has reached a steady 
state. The time independent value of c, in the frame of 
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FIG. 3. Net force on the last boat vs. its distance from the 
next forward boat. A positive force will decrease the distance. 



the boats is given by eqlS) 
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dx 
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The solutions to this equation for one boat and two 
boats are shown in fig 12b. Note the higher concentration 
around the trailing peak. This feature remains even if the 
boats are placed close enough to be in contact. By the 
mechanism described above, the two boats will approach 
each other and travel in a close, stable configuration. If 
additional boats are added to the system, they will be 
gathered into the cluster through the same mechanism. 
This leads to the question of whether the cluster remains 
stable as the number of boats increases. 

Using the constant velocity assumption, it is possible 
to construct a stable cluster with many boats. We can 
then set the location of the rearmost boat and calculate 
the net force on it as a function of distance from the next 
forward boat. This allows us to determine the stability 
of the cluster. 

Under certain conditions, the resulting net force shows 
a stable equilibrium point at a relatively close distance 
from the next boat (see figE]). The region of attraction, 
shown as an interval with positive net force, is long com- 
pared to the boat length. This makes the equilibrium 
stable against small perturbations and ensures that boats 
could be drawn in from an appreciable distance. If the 
net force remains positive at zero spacing, the boats will 
remain in contact as they travel. It is important to re- 
member that this mechanism only applies while the aver- 
age camphor concentration around each boat is less than 
Ct . If c is higher than this value, the opposite effect will 
slow down the trailing boat and effectively spread the 
boats apart. The periodicity of the system also places 
a limit on the number of boats that can form a stable 
cluster as explained below. 

Comparing these results with the observations of ants 
reveals many similarities. As described by John et al. 0, 
a line of ants may form platoons. These stable, high den- 
sity groups of ants move at a constant speed. Individual 
ants following behind the platoon sense its pheromone 
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FIG. 4. a) b changes abruptly at particular values of k and 
a. b) There is a maximum value of A'', dependent on various 
parameters, above which clusters do not form. The other 
dimensionless parameters used here are fi = 0.24, F — 528, 
R = 75.8, ro = 0.25. 



trail and may move at a higher velocity. This allows the 
individual ants to catch up with and join the platoon. 
Clusters of camphor boats are also stable, high density 
groups that travel at constant speed. The trail of cam- 
phor left by the cluster causes following individuals to 
travel at higher speed and catch up with the cluster. Al- 
though camphor boats are not sentient and are driven by 
much simpler principles, these behaviors are very similar. 



IV. NUMERICAL RESULTS 

It is important to provide a quantitative way to discern 
clustering from other kinds of behavior. Considering the 
structure and stability of clusters, we propose the follow- 
ing quantity as an order parameter. 



max 

i=l..N 



R- NL 



(9) 



For boat number i, dfi and dti are the distances from 
the next boat in front and in back respectively, and the 
normalizing factor, R — NL, is the total empty space in 
the route. The brackets represent a time average. For a 
system containing one cluster with all boats in contact, 
6=1. For ideal homogeneous flow, 5 = 0. For jammed 
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flow, b will approach zero when averaging over long time 
intervals. In practice, particularly in experiment, b may 
deviate slightly from these values, but should generally 
be much larger for the clustered mode. 

By testing the dependence of b on the model parame- 
ters, we found that the most important dependence was 
on the parameters related to the camphor level. These 
include the evaporation constant, k, and the supply pa- 
rameters, a and tq . The relation to k and a can be seen 
in Fig|3h.. Since the effect of vq is very similar to that of 
a, and its values are very limited in practice, it was not 
examined separately. 

There was also a dependence on the number of boats. 
One reason is that adding more boats will increase the 
total camphor being added to the water. Also, if the 
size of the cluster becomes comparable to the total route 
length, the boats at the back of the cluster will influence 
the boats at the front. In simulation this effect disrupted 
the stability and destroyed the cluster. Considering these 
two effects, there appears to be a maximum number of 
boats above which clusters do not form as shown in fig HI 
This maximum number depends on several factors, but 
our tests suggest that it is typically less than the critical 
number required for the jammed mode to appear. Thus, 
for the parameter space tested, clustering and jamming 
did not coexist. However, it is possible that there exists 
some set of parameters for which clustering and jamming 
can both occur. 

Considering the very abrupt changes in b, it is worth- 
while to check for hysteresis. The results shown in figE] 
were obtained by gradually changing either the number 
of boats N or the supply rate a. The system was given 
time to approach a steady state between each measure- 
ment. The figure shows clear hysteresis with respect to 
the number of boats and the supply rate for this set of 
parameters, but if k is increased, the hysteresis with re- 
spect to a vanishes. This may suggest a transition from 
a subcritical to a supercritical bifurcation. 



V. CONCLUSION 

To better understand clustering behavior seen in ex- 
periment, we examined a mathematical model of an en- 
semble of camphor boats in a periodic, one-dimensional 
system. We found that this behavior was caused by the 
relationship between surface tension and camphor con- 
centration on the surface of water. For low camphor 
concentrations, boats may experience a stronger driving 
force by traveling close behind other boats. We then de- 



fined an order parameter to quantify the clustering be- 
havior and applied it to the results of numerical simula- 
tions to examine the dependence on system parameters. 
Our next goal is to investigate this phenomenon exper- 
imentally to test the validity of the theoretical analysis 
presented here. Clustering behavior can be reproduced 
in experiment, but the dependence on experimental con- 
ditions has not yet been determined. 
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FIG. 5. Hysteresis is seen by gradually increasing and de- 
creasing N (k = 0.73) and a {k = 0.36). For k = 0.73, the 
hysteresis with respect to a vanishes. Other parameters as in 
FigE 

In analyzing this model of collective motion, we hope 
to have added some insight into clustering or swarming 
behavior. The striking similarity to the behavior of ants 
seems to support the idea that collective motion can be 
categorized into a set of universal classes. Although the 
mechanism presented is specific to the motion of cam- 
phor boats, it represents yet another way in which nature 
forms organized structures from relatively simple pieces. 
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